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ABSTRACT 

We present results from modeling of quasi-simultaneous broad band (radio through X- 
ray) observations of the galactic stellar black hole (BH) transient X-ray binary (XRB) systems 
XTE Jl 1 18+480 and GX 339-4 using an irradiated disc + compact jet model. In addition to 
quantifying the physical properties of the jet, we have developed a new irradiated disc model 
which also constrains the geometry and temperature of the outer accretion disc by assuming 
a disc heated by viscous energy release and X-ray irradiation from the inner regions. For the 
source XTE Jl 1 18-1-480, which has better spectral coverage of the two in optical and near-IR 
(OIR) wavelengths, we show that the entire broad band continuum can be well described by 
an outflow-dominated model + an irradiated disc. The best-fit radius of the outer edge of the 
disc is consistent with the Roche lobe geometry of the system, and the temperature of the 
outer edge of the accretion disc is similar to those found for other XRBs. Irradiation of the 
disc by the jet is found to be negligible for this source. For GX 339-4, the entire continuum 
is well described by the jet-dominated model only, with no disc component required. For the 
two XRBs, which have very different physical and orbital parameters and were in different 
accretion states during the observations, the sizes of the jet base are similar and both seem to 
prefer a high fraction of non-thermal electrons in the acceleration/shock region and a mag- 
netically dominated plasma in the jet. These results, along with recent similar results from 
modeling other galactic XRBs and AGNs, may suggest an inherent unity in diversity in the 
geometric and radiative properties of compact jets from accreting black holes. 

Key words: black hole physics — radiation mechanisms: non-thermal and thermal — X-rays: 
binaries — stars: winds, outflows- accretion, accretion discs- X-rays: individual: (GX 339-4, 
XTE Jl 118-1-480) 



1 INTRODUCTION 

Soft X-ray transient (SXT) binary systems that occasionally go into 
episodes of high activity or outbursts are some of the best labo- 
ratories to study accretion flows near compact (GM/R ~ c^) ob- 
jects. Coordinated multi-wavelength campaigns have been carried 
out in recent years to observe SXTs simultaneously over a broad 
frequency range of nine decades or more from radio to X-rays 
(see, e.g.. iFendeil [20 06. for a review). These campaigns have pro- 
vided extremely valuable insights about the accretion in/outflows 
and helped to constrain the theoretical models predicting accretion 
geometry and physical emission processes. Outbursts of SXTs are 
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usually expl ained in terms of a disc instability model (see, e.g., 
lLasotall200ll . for a review), postulating a dramatic increase in mass 
accretion via an accretion disc towards the central compact object. 
The luminosity of a black hole transient (BHT) can vary by as much 
as eight orders of magnitude between quiescence and peak outburst, 
within timescales of weeks to months. Such systems are also ideal 
for studying the formation of jets, which are a by-product of accre- 
tion, yet poorly understood so far. 

Matter spiralling into the gravitational potential well of the 
compact object forms an accretion disc, and radiates via dissipation 
of the accumulated gravitational potential energy. The observed 
UV and soft X-ray spectra of these sources, during epochs of high 
mass accretion rate, are well modeled by viscous dissipation from 
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the inner regio ns of a geometrically thin, o ptically thick "Shakura- 
Sunyaev disc" jShakura & Sunyaevll 19731) . 

However in the OIR, the viscously dissipating disc model 
sometimes significantly underpredicts the observed spectral 
energy distribution (SED), and the excess is attributed to irra- 
diation of the outer disc by the intense X-ray radiation from 
the inner region close to the compact object. Th is irradiation 
mech a nism has been we ll explored t heoret ically (Cunningham 
19761; IVrtilek et al.l 119901 : iKing et al.l 1 19961 : iDubus et al.. .1999. , 
200 ih and thei r eff^ect obse r ved in optical and near-IR wave- 
lengths (e.g. Vrtil ek et al. ' '1990'; 'van Paradijs & McClintock! 
19941: iMcClintock et al. 1995; Hynes et al. 2002; Hvnes 2005 



Russell et alJl2006l : iMaitra & Bailvnll2008i : loierlinski et alJI2009h . 



The most prominent, observable effect of irradiation is to heat up 
the outer regions of the accretion disc to temperatures higher than 
what would be possible solely by viscous dissipation of a standard 
Shakura-Sunyaev disc. This leads to a deviation in the shape of the 
SED of the accretion disc (from that of a standard disc), creating 
a 'bump' in the optical through infra-red (OIR) wavelengths. 
As the quality of the OIR data improves, the deviation (i.e. the 
bump) from Shakura-Sunyaev spectrum becomes more and more 
apparent, and the need to incorporate irradiation in the model 
becomes important. 

In Roche lobe overflowing BHTs where the accretor is 10 
times or more heavier than the donor. X-ray irradiation of the 
donor is small. Furthermore irradiation causes the out er disc to 
flare, shielding most of the donor star from irradiation dMilgroml 
Il978l) . While most models assume an azimuthally symmetric accre- 
tion d isc, warping of the accretion disc (see, e.g.. lOgilvie & DubusI 
I2OOIL and references therein) can introduce additional uncertaini- 
ties in shielding the donor from irradiation. 

Not all the mass accreted inwards falls onto the compact ob- 
ject. Outflows, in the form of steady compact jets, or transient, bal- 
listic ejecta are often seen in compact binary systems (see, e.g., 
[Fender 2006, for a recent review). Steady, compact jets have been 
observed from an increasing number of compact binary systems 
when they are in an X-ray hard state (see below). The hallmark 
signature of such a jet is its fla t to slightly inverted spectrum which, 
following the classic work bv lBlandford & Konigll ( [l979h . is usu- 
ally attributed to a superposition of self-absorbed synchrotron emis- 
sion from segments of a collimated jet. For stellar BHTs the flat 
spectrum of the jet is ex pected to extend fro m radio frequencies 
up to IR (~ 10" '" Hz: iMarkoff et al.ll2003h . beyond which the 
jet sy nchrotron becomes optically thin (see, e.g.. lCorbel & Fendej 
I2OO2) . At X-ray frequencies, inverse Compton (IC) upscattering in 
the jet base of synchrotron as well as external photons from the ac- 
cretion disc can become com parable and m ay even dominate the 
spectrum. As shown in Mark off et al.l ( |2005|) . synchrotron and IC 
at the jet base can result in spectra very si milar to those produced 
by co mpact Comptonizing coronal models ( |Poutanerj[T998l : ICopDil 
1 19991) . 

During the course of an outburst, the X-ray spectral and 
temporal features show an enormous amount of diversity. How- 
ever extensive observations and their analyses over the past few 
decades have enabled the classification of spectral and temporal 
feat ures into relatively few groups known as X-ray states ( see, 
e.g.. iMcClintock & Remill"^ l2006l : iHoman & BeUonill2005l, for 
extensive definitions and review of X-ray states). In particular, 
two canonical states are often seen: a "soft" state characterized by 
strong thermal emission from the accretion disc and low variability 
in the light curve, and a "hard" state characterized by non-thermal 
power law emission, high variability and probably a collimated out- 



flow. Besides these canonical states, a host of intermediate and ex- 
treme states with varying ratios of thermal and n on-thermal radia- 
tion also are seen. Recent OIR observations (e.g. lMaitra & Bailvnl 
2008; Russell et al. 2008) covering complete outbursts of many 
XRBs have shown that while irradiation usually plays a dominant 
role in the evolution of the OIR spectrum during the entire out- 
bursting phase, the contribution from a non-thermal source, e.g. a 
jet, can become dominant when XRBs (most notably black holes) 
are in hard state. Thus any model which aims to understand disc-jet 
coupling must include irradiation. 

We have developed a new spectral model which calculates 
continuum emission from an irradiated accretion disc, and in con- 
junction with a modified ve rsion of the compact jet model devel- 
oped by Ma rkoff et alj 12005.) , allows us to constrain the geometri- 
cal and radiative properties not only of the jet, but simultaneously 
those of the accretion disc. The thermal photons from the irradiated 
disc are also included in the photon field of the jet for Compton 
scattering, although due to Doppler redshifting in the jet frame, this 
external Compton emission is usually much smaller than the syn- 
chrotron self Compton (SSC). The new model is fully integr able 
with the standard spectral analysis software ISIS jHouckl2002l) and 
parallelizable for using in large cluster computing envirormients 
(also see 33.3b . 

We have used this new model to analyze broad band 
quasi-simultaneous observations of two galactic BHT sys- 
tems XTE Jl 118+480 and GX 339-4. Simultaneous or quasi- 
simultaneous observations of both sources, during hard X-ray state 
of their outbursts, have revealed broad band continuum emission 
rang ing from radio to X-rays, strongly sug gesting th e presence of 
jets (iMarkoff et alj|200ll l2003l : iHoman et al. 2005; Markoff et all 
2005). By studying short time-scale variability in near-IR wave- 
lengths from XTE JI 118+480, .Hynes et al . (2006) also suggested 
a non-thermal origin of the near-IR radiation, with a power law in- 
dex typical of optically thin synchrotron emission from a jet. 

The known physical properties of the two sources 
XTE JI 118+480 and GX 339-4 have been briefly summa- 
rized in ^ Since analysis of multi-wavelength data requires 
careful calibration of data from various ground-based as well as 
space-borne instruments, we have also presented the data reduction 
and calibration procedures in ^ A brief outline of the model, with 
somewhat more emphasis on the newly introduced irradiated disc 
component, is given in ^ along with a description of the fitting 
procedure. The results of our broad band SED modeling suggest 
that despite the differences in size, mass, orbital parameters or 
accretion state, some of the fundamental physical parameters 
characterizing the geometrical and radiative properties of the 
jet are similar for both sources. This similarity may be a global 
property of jets formed near compact objects and is discussed in 

m 



2 SOURCES, OBSERVATIONS & DATA REDUCTION 

XTE Jl 1 18+480 and GX 339-4 are two prototype galac- 
tic BHTs with markedly different properties. Both binary sys- 
tems contain hi gh mass compact objects, 8.5 ± O.6M0 for 
XTE Jl 1 18+480 loelino et al.ll2006l) and > 6 Mq for GX 339-4 
dMunoz-Darias et alj|2008l) . thus strongly suggesting for the com- 
pact objects to be black holes in both cases. 

XTE JI 1 18+480 lies 62°. 3 north of galactic plane, i n the halo. 
It has an orbital period of 4.04 hours l lGelino et alj|200^ . The dis- 
tance and interstellar column density to XTE Jl 1 18+480 are rea- 
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sonably well known from OIR photometry dur ing quiescence to b e 
1.72±0.1 kpc and 1.3x lO^" cm'^ respectively jOelino et alj2006h . 
We adopted a black hole mass of 8.5 M© for this sour ce. The orbita l 
inclination of XTE Jill 8+480, as measured by G elino et al.l j2006l) 
is 68 ± 2 degrees. 

GX 339-4 lies near the galactic plane, only 4°. 3 south of 
the galactic equat or and has an orbital period of 42.14 hours 
( lHvnesetalJl2003h . We adopted a black hole mass of 7 Mq for 
GX 339-4 in this work. The secondary star in the GX 339-4 sys- 
tem has eluded detection so far, thus making estimates of masses, 
distance or line of sight extinction difficult. Here we have used a 
distance of 6 kpc and A'h = 6x 10-' cm~^ for GX 339-4 whic h are 
consistent with the limits given by iHvnes et all ilOOA |2004|) and 
iMuiioz-Darias et alj ( 12008) . 



2.1 OIR and radio data 

XTE Jl 118+480 went through an outburst during January- 
February 2005. A discussion of the OIR light curve morphology 
of the 2005 OIR outburst of XT E Jl 118+480 is p resented by 
IZurita et al.l \200dt) . It was noted bv lZurita et"al] ( l2006h that the 5- 
12 keV/ 3 -5 keV hardness ratio from the All Sky Monitor (ASM; 
iLevine era l. 1996) data onboard the Rossi X-Ray Timing Explorer 
(RXTE; lBradt et al.lll993) satellite suggested that the source was 
in an X-ray hard state during the entire outburst of 2005. For this 
source we used quasi-simultaneous optical (from the Liverpool 
telescope), near-IR (from UKIRT) and 15 GHz radio data (from 
the Ryle telescope). The full set of multi-wavelength observations 
obtained during this outburst of XTE Jl 1 18+480 will be presented 
in Brocksopp et al. (in prep). 

XTE Jl 118+480 went through an outburst in 2000 as well. 
During this outburst the source was extensively observed by many 



ground based telescopes as well as space-borne missions (see e.g . 
Hvnes et al.l I2OO0I : iMcChntock et alj I2OOIL iFrontera et al] I2OOII : 



Chatvetal.H2003l) . In particular, the source was observed quasi- 



simultaneously around 2000 April 18th using the Ryle radio tele- 
scope, UKIRT, HST, EUVE, Chandra and RXTE giving an un- 
precedented broadband spectral co verage. This excellent data set 
has been used in severa l works (e.g. lMarkoff et alj200ltlEsin et al.l 
l200lMYuan et alj2005l) to infer properties of the source. We revisit 
this data set to constrain the geometry and radiative properties of 
both the jet and the disc using our new disc+jet model. 

We used th e OIR and radio SED of GX 339-4 reported 
by iHoman et al. I l l2005l) . which were obtained during a multi- 
wavelength campaign to observe its outburst of 2002. For all data 
sets we used the photometric OIR magnitudes and de-re ddene d 
them using the interstellar extinction law of ICardelli et alj ( Il989h . 
In the absence of definitive measurement of the ratio of total to se- 
lective extinction (i?v[= Av/E{B — V)]) in the direction of either of 
the sources, we have taken Ry = 3.1, the widely used value mea- 
sured by Rieke & Lebofsk^ j 19851) towards the galactic center, and 
standardly used for these sources. Thus it must be kept in mind that 
if the "true" value of Ry is significantly different from 3.1 (given 
neither of the sources are very close to the galactic center), it might 
introduce additional uncertainity in reddening correction. 



2.2 X-ray data 

For both so urces, we used data from the Proportional Counter Ar- 
ray (PCA; IJahoda et al]|20()6l) and the High Energy X-Ray Tim- 
ing Experiment (HEXTE; Rothschild et al.l 1 19981) onboard the 



RXTE satellite. The X-ray data were downloaded from the NASA 
HEASARC's public archivtfl and reduced using an in-house script 
following the standard reduction procedures outlined in RXTE cook 
fcooy(0 using HEASOFT software (v6.1.2 Arnaud 1996). To avoid 
any terrestrial contamination, we considered data only for eleva- 
tion angles of 10° or higher. Also any data with pointing offset more 
than 0°.02, or within 30 minutes of SAA passage, or trapped elec- 
tron contamination ratio more than 0.1 were rejected. Since both 
the sources were bright (PCU2 count rate > 40 counts/s/PCU), 
the model pcajbkgd-cmbrightvle_eMv20051128.mdl was used for 
background subtraction. Data from all the layers of the active PCUs 
during the observations were extracted to create the spectra. Sys- 
tematic uncertainty of 0.5% was added in quadrature to the PCA 
data to account for uncertainties in the calibration. For the HEXTE 
instrument, data from both clusters were added. Both PCA and 
HEXTE data were binned to achieve a minimum S /A' of 4.5. For 
the PCA we used the energy range of 3-22 keV and for HEXTE 
20-200 keV. The normalization of the OIR and radio data was tied 
to that of PCA. HEXTE normalization was allowed to vary. 



3 MODEL, ANALYSIS AND RESULTS 

The continuum jet model used in this work is based on 
iMarkoff et al.l j2005l) , (henceforth referred to as the agnjet 
model). Given the importance of irradiation in the observed spectra 
of XRBs, we have modified the model to include additional physics 
(see 33.2l for details) that allow us to calculate the effect of irradia- 
tion on the outer disc due to high energy photons from (1) the jet, 
and (2) the inner regions close to the compact object, impinging on 
the disc. We also compute inverse Compton scattering of thermal 
disc photons by the electrons in the jet plasma, although their effect 
on the total spectrum is negligible, except near the jet base, due to 
Doppler redshifting of the disc photons in the jet frame. 



3.1 Jet parameters 

The details of the physics of the jet model as well as a full de- 
scr iption of its main in put parameters is given in the appendix 
of lMarkoffet"ai] j2005l) . Therefore we only give a brief sum- 
mary h ere, and outline the modifications made to the lMarkoff et al.l 
( |2005|) model: 

The main parameters that determine the properties of the jet 
are the input jet power (A',), electron temperature of the relativistic 
thermal plasma entering at the jet base (T^), the ratio of magnetic to 
particle energy density (a.k.a. the equipartition factor fc), physical 
dimensions of the jet base (assumed to be cylindrical with radius 
To and height h^), and the location of the point on the jet {Zacc) be- 
yond which the a significant fraction of the leptons are accelerated 
to a power law energy distribution. Nj, parameterized in terms of 
the Eddington luminosity, determines the power initially input into 
the particles and magnetic field at the base of the jets. In the ab- 
sence of a full understanding of the mechanism for energizing the 
jets, Nj plays a similar role t o the compac tness parameter in ther- 
mal Comptonization models ( ICoppilll999l) . A typical value of the 
number density of leptons (n) at the base of the jet (derived from 
fits) is n = N'^L^iiKAy^p^m^c^nr^) ~ lO'''"'' cm"'. Here Le^j is the 
Eddington luminosity of the source, ji^ ~ 0.4 is the sound speed at 
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the base (see below), - (1 -/Jj)"''^, is the proton rest mass, 
and To ~ 10' cm is the radius of the jet base, implying an optical 
depth T = najra ~ 0.001 - 0.01. Since the optical depth is small, 
the probability of multiple Compton scattering is very small and 
we only consider single scattering. For the non-thermal power law 
electrons, the lower cutoff in the Lorentz factor is assumed to be 
equal to peak of the relativistic Maxwell-Juttner distribution, i.e. 
7e;„„„ = l.l'ikTJim.c^). 

Followin g the prescription of iFalcke & Biermamil ( Il995h : 
iFalckd ( Il996h for maximally efficient jets, we assume that the bulk 
speed of the plasma (with adiabatic index F = 4/3) at the base of 
the jet is the sound speed given by = V(r - 1)/(F + 1) ~ 0.4. We 
do not model the particle acceleration in the jet, but just start at the 
sonic point. Each jet then accelerates longitudinally along its axis 
due to pressure gradient and expands laterally with its initial sound 
speed. The bulk velocity profile, magnetic field strength, electron 
density and electron Lorentz factor along the jet are calculated by 

solving the adiabatic, relativistic Euler equation. 

In the previous versions of the jet model (e.g. in lMarkoff et al.l 

l2005h we used the ratio of scattering mean free path to gyroradius 
of the particles in the plasma through which the shock moves (/,r), 
as a model parameter. However given the current lack of physical 
understanding of the actual shock propagation mechanism, in this 
paper we fixed the shock speed in the plasma (/?,;,) to 0.6 and used 
the quantity e^^ — /(ratio of scattering mean free path to gyro- 
radius) as a model parameter instead of The parameter e^c can 
be physically interpreted to be a param etrization for the sho ck ac- 
celeration rate (fa„- <^ f^,'; see eq. 1 of lMarkoiret"ai]|200lh . The 
dominant cooling processes within the jet are synchrotron and syn- 
chrotron self-Compton radiation. If the inner edge of the accretion 
disc is sufficiently close to the base of the jet, then the soft photons 
from the accretion disc, acting as seed photons for external Comp- 
tonization, can also cool the particles in the jet. The jet continuum 
as observed by an observer on Earth is a superposition of spectra 
from individual segments along the jet axis (taken to be perpendic- 
ular to the accretion disc plane), calculated by solving the radiative 
transfer equation, after taking into account relativistic beaming ef- 
fects. 



3.2 The irradiated disc 

Since the coupling between disc and jet is not well understood, 
we assume an independent classical thermal viscous "Shakura- 
Sunyaev" accretion disc, with a radial temperature profile T(R) oc 
i?"''"*, and modeled by the parameters Ti,,, rj,,, the temperature and 
radius at the inner edge of the disc. If the outer parts of the disc are 
irradiated, then as discussed below, it is further possible to constrain 
the temperature and radius at the outer edge of the disc. Given the 
mounting evidence for the existence of relativistic outflows from 
black holes, we also explore the influence of the jet acting as a 
source of irradiation heating of the accretion disc. 



3.2.1 Can jets influence the energetics of the accretion disc ? 

We computed the radiative jet flux incident on the accretion disc as 
a function of radial distance on the disc plane, taking into account 
special relativistic beaming effects. In order to obtain an upper limit 
on the jet-induced heating of the accretion disc, we assumed that 
the entire incident radiative flux from the jet fv(R) is thermalized 
upon being absorbed by the disc, thus heating the disc locally at 
radius R ala. rate Hj^tiR), given by 



H,AR)= f MR)dv (1) 
Jo 

We calculated Hj^i(R) at ten logarithmically spaced radii span- 
ning from Win to /?out- The ratio of /fjet(W) to the local viscous heat- 
ing (//„,sc(W)) is shown in Fig.Q At large radii (R > 100 Rg), Hj^iR) 
has a ~ r^^ * radial profile, which is somewhat steeper than that ex- 
pected from irradiation by a static corona or the inner accretion 
disc (for which the irradiatio n heating falls as ~ r"'^ '° "^"i see, 
e.g., O!22]and|Hsnei|2003)). However the local viscous heating 
rate drops outward even faster (Hyisc(R) R^), thus causing the ra- 
tio Hja{R)/H^isdR) to slowly increase for increasing R. As the disc 
radius becomes small and approaches that of the jet base, the solid 
angle subtended by the jet increases rapidly causing the jet heating 
to saturate. However the viscous heating continues to rise at smaller 
radii and causes /fji.,(/?)//fvisc(^) to drop at small radii. However, as 
evident from Fig.[T] even the maximum possible jet-induced heat- 
ing of the disc is about seven orders of magnitude smaller than the 
viscous heating. This confirms that relativistic beaming of the jet 
causes a very small fraction of the flux emitted by the jet to fall 
back on the disc, and the resulting influence of the jet on the ener- 
getics of the disc is negligible. 

Light bending effects near the black hole can enhance the flux 
incident on the disc, but by not more than a few tens of percent of 
the numbers calculated above (i.e. by assuming Minkowski metric 
instead of a more realistic Kerr or Schwarzschild metric). 



3.2.2 An irradiation source near the compact object 

Analy sis of OIR fight curves of XTE II 118-1-480 by IZurita et ahl 
( l2006h suggests there might have been significant contribution of 
thermal flux in the optical. The near-IR flux however is reported 
to be significantly non-thermal in nature from studies of rapid 
time variability in the light curve, a s well as nearly flat IR SED 
jZurita et alj2006l ; lHvnes et al.l2006t) . As discussed in the previous 
paragraph as well as shown in Fig.[T] the jet's contribution in heat- 
ing the disc is negligible; we have therefore added an irradiation 
heating term due to a somewhat more static source of irradiating 
X-rays near the black hole. The physical origin of this irradiating 
source could be the inner accretion disc near the black hole. Since 
we assume a radial temperature profile and do not solve the local 
disc structure (which in itself is a detailed MHD problem, see, e.g., 
.Hawley et al. 2007, and references therein), our model cannot dis- 
cern the source of the irradiating X-rays. 

The temperature due to irradiation heating (rinad) is usually 
assumed to have a power law radial dependence of the form R^". 
Depending on the initial assumptions about disc structure and ge- 
ometry of the irradiating s ource, theoretical models predict n to 
be in the range o f 0.4-0.5 jVrtilek et al]|l990l : iKing et al.l [19971: 
iDubus etalj|l999l) . Given the uncertainities in de-reddening and 
scarcity of data points in the OIR, the exact choice of n is not vi- 
tal for a global understanding of the broad-band SED. We used 
w = 3/7 (vertically isothermal disc with disc heig ht h oc r^^'), 
mainly for consistency with earlier work (e.g. by IVrtilek et al.l 
ll99(]|jHvnes et alj2002l : lHvnesll2005h . 

We assume that the total disc heating is a sum of local viscous 
heating, and heating due to a static source of irradiation near the 
black hole. The effective temperature of the disc (r^ff) is therefore 
related to the viscous temperature (ryisc) and the irradiation tem- 
perature (riirad) as follows: 

Tts(R) = Tt,JR) + TljR), (2) 
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where T,,, = T,,(R/R,J-^'* and T,^,^ = TURI RoutY^" with 
T'in/outi ^in/out as free parameters. Other parameters, viz. masses, in- 
chnation and distance to the binary system, are taken from values 
published elsewhere and kept fixed during the process of fitting. 
The radius /fimd where irradiation heating becomes equal to vis- 
cous heating can be estimated from 

Rmnd = \('rmlToud R-ali'' IR\n'''\ ' 

Since viscous dissipation dominates for R < i?i,iad, the disc height 
at j ?irrad is estimated from Shakura-Sunyaev's a -disc solution (see 
e.g. lShakura & Sunvaevlll973l : lFrank et aPl 20021) as 



H (i?,,,,d) = 1.7 X 10« ff-'/'o M]l^-\MIM^r^' i?r,;:d;io t 



-3/8 r,9/8 



f3/5 



(4) 



where we have Mj^ = M/CIO'* g/s), R\n-iA:\o = 
^mad/CIO'" cm), / = [1 - 2GM/(c-2i?,,,ad)]"'' and chosen a = 0.05. 
Similarly, since H oc ij''' when R > i?i,iad, the disc height at i^out is 
given by 



H (i?oiit) - ^(^iirad) X (^out/^inad) 



9/7 



(5) 



The solid angle subtended by the outer, irradiation dominated disc 
as seen from the compact object is computed from equations |4] and 

s 

The agnjet model computes the broadband continuum from 
the jet and the irradiated disc but not line emission or reflection 
features. Therefore we used an additional Gaussian line profile to 
fit the iron Ka e mission complex near 6.5 keV a nd the convolution 
model reflect jMagdziarz & Zdziarskj||T99^ for the excess in 
the range of 10-30 keV, both features being usually attributed to 
reflection by the relatively cold accretion disc. The column density, 
taken from elsewhere published values, was fixed during the fits. 
For the reflect model, which accounts for reflection from neutral 
material in the accretion disc, we used solar abundance and tied the 
inclination to that of the disc (and also the jet axis). The energy of 
the iron line in the Gaussian model was allowed to vary between 
6-7 keV and its width was fixed to 0.5 keV. The spectral coverage, 
along with column density (A^h) towards the sources, distances, and 
masses of the compact objects are given in Table[T] 



3.3 Model fits 



Spect ral fittings were performed using ISIS (v 1.4.9-4; iHouckl 
l2002h . Evaluation of confidence intervals for any physical model 
is computationally expensive, and this is true for agnjet model 
as well. Therefore we have used the parallelization technique de- 
scribed in lNoble et al. Ill200d) and distributed the task of computing 
confidence interval^ for / free parameters over //2 or (/ + l)/2 
nodes (depending on whether / is even or odd; each nod e is an Intel 
Xeon 3.4GHz processor) of a Parallel Virtual Machine dOeist et al] 
Il994h running on the LISA cluster in Almere, Netherlands. This re- 
duced the runtime to under 24 hours, a greater than 75% speedup; 
it also allowed us to discern finer features in parameter space by 
increasing the tolerance resolution by a factor of 500, while keep- 
ing the overall runtime to approximately half the serial runtime of 
96 hours. While these are welcome improvements, reducing the 
model runtime even further would enable us to analyze new data 
more quickly. For this reason we are experimenting with ways of 
utilizing more than one processor per parameter, such as partition- 
ing parameter space more finely, and exploring the use of OpenMP 

^ Using our own custom versions of the cljnaster and cl.slave scripts 
described at |http://space.mit.edu/cxc/isis/single-param_conf .limits .html| 



jPagum & MenonI 1 1998b to parallelize internal loops within the 
model source code. 

Using the distributed parallel computing technique described 
above, we determined the best-fit parameter values as well as their 
confidence intervals. Our confidence intervals correspond to l^x^ = 
2.71 for a given parameter (which for normal distribution would 
imply a 90% single parameter confidence limit). 

3.3.1 XTE Jl 118+480 during 2005 outburst 

Of the three data sets presented in this work, the data from the 2005 
outburst of XTE Jl 1 18+480 was the faintest and had a steep power 
law slope in the X-rays. The flux in the HEXTE bands was ex- 
tremely low for this data and while regrouping the data we found 
that no flux with S/N > 4.5 could be detected above 60 keV. Since 
the single point in the radio frequency does not constrain the length 
of the jet, we fixed it to 10" ' cm which is a lower limit on the ra- 
diative length of the jet assuming that the radio data point is still 
on the optically thick synchrotron regime with a flat-to-slightly- 
inverted spectral index. While the fits require a high fraction of 
non-thermal electrons in the post-shock jet (> 0.5), they are not 
very sensitive to the exact value of this fraction, and we fixed its 
value to 0.75. The data cannot be described by a non-thermal jet 
continuum alone and an excess in the OIR data requires additional 
flux which we modeled using the irradiated accretion disc model 
described in 93.21 This thermal excess in the optical h as also been 
suggested from the OIR cove rage of the outburst by IZurita et al.l 
( 200^ and lHvnes et al] ( |2006|) . 

The X-ray data (3-60 keV) can be well fit (x^/v = 88.8/55) 
only with a power law, however systematic residuals in the 5-7 keV 
range suggest that a weak iron line may be present. Adding a Gaus- 
sian line to the power law improves the fits tox^/y - 79.5/52, with 
the line center near 6.6 keV and width of 0.1 keV. Given the weak 
strength of the lines, we fixed the line energy as well as the width 
to these values, and did not use any reflection model modeling the 
broad band SED of XTE Jl 1 18+480. 

Fits performed by allowing rjn and 7,^ of the accretion disc 
to vary freely results in good fits with quite small values of both 
Tin and Tin. There are uncertainities in our understanding of the 
physical properties of the innermost regions of the disc close 
to ISCO, dereddening of interstellar extinction, and also uncer- 
tainities in instrumental response and calibration. These concerns 
make detection of a cold disc with small r^, Tjn and its physi- 
cal origin o pen to debate (Miller et al."2006a''b';'Rvkofi" et al. 2007|; 
Ramadevi & Seetha ,2007 ; Gierlihski et al. ,2008, ; D'Angelo et aU 



20081; iReis etai]|2009t) . Another issue with the XTE Jl 1 18+480 



data is the jet/disc inclination. While we usually fix the inclina- 
tion to a value published elsewhere, the extremely flat radio to IR 
spectrum of XTE Jl 118+480 requires a much smaller inclination 
than the published value of ~ 70° (see, e.g., Gelino et al. |2006l ; 



IChatv et al.ll2003l ; iMcClintock etai]|200ll ; IZurita etalj|2002h . Al 
lowing the inclination to vary for this source results in jet inclina- 
tion of ~ 25°. 

Therefore we tried two different models for this data set of 
XTE Jl 118+480. In the first model (Model 1; second column in 
Table |2j, apart from other regular parameters described above, the 
jet/disc inclination as well as r^^ and T]„ were allowed to vary freely. 
In the second, somewhat more conservative model (Model 2; third 
column in Table|2j, the jet inclination was fixed to 30°, ri„, Tjn were 
also fixed to 30 Rg and 0.1 keV respectively. The best-fit values 
for the second model also gives similar jet parameters as the first 
model, although with a slightly higher jet power and steeper elec- 
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tron distribution index. Also the second model favors a somewhat 
smaller and hotter outer disc edge. 



3.3.2 XTE Jl 118+480 during 2000 outburst 

We used the broadband SED data of XTE Jl 118+480 obta ined 
around 2000 April 18, published in iMcClintock et alj j200lh and 
iHvnes et al. While analyzing this data set a broad dip in the 

combined EUVE and Chandra spectra was noted between 0. 15-2.5 
keV (see e.g. the residuals in Fig.[3]i which cannot be modeled using 
the standardly used column density for this source ( 1 .3x 10^^ cm"^). 
This feature was been noted by previous works as well, and is 
attributed to me tal absorption in a partially ionized gas, (see e.g. 
lEsin et al]|200ll) . We have therefore excluded the energy range be- 
tween 0.15-2.5 keV from our broadband fitting. Results of fitting 
indicate that during this observation emission from the jet dom- 
inates the radio and IR regions. Contribution from the accretion 
di sc starts becoming d ominant in the optical. However, as discussed 
in lHvnes et al.l ( I200Q ). the EUVE fluxes are extremely sensitive to 
the extinction law, and the slope of the dereddened EUVE SED 
can range between +2 > a > -4 (F,, ~ y") for the permissi- 
ble r ange of Nh . Here we have used the EUVE fluxes reported by 
( lH7nes e t al. 2000; McClintock et al. 2001) and a column depth of 
Nh = 1.3 X 10"" cm"^, which results in a very steep EUVE slope. 
This cutoff in the EUVE data would imply a cold, truncated inner 
disc with with ri„ ~ 340 R^ and Tin ~ 2.9 x 10^ K. Given the low 
temperature of the inner disc and the consequent absence of a ir- 
radiating soft X-ray photons, it is not unexpected that we do not 
find any signature of irradiation in this case (i.e. T„^, is not con- 
strained). As in the case of the 2005 data of this source, the X-ray 
emission in this model is entirely dominated by optically thin syn- 
chrotron emission from the jet. While the fits require a small nozzle 
(hg/rQ ~ 0.2 - 0.6), they are not very sensitive to the exact value 
of this parameter, and we fixed its value to 0.4. We note that this 
value is abot a factor of three smaller than what is found typically 
in other sources, or even during the 2005 outburst of this source. 



3.3.3 GX 339-4 during 2002 outburst 

The X-ray data of GX 339-4 cannot be well fit by a power law 
alone. Using simple power law+ Gaussian model gives unaccept- 
able fit to the data with reduced of 5.9. Convolving a power- 
law-l-gaussian with a reflection model give better fits, but still with 
rather a large reduced of 2.3. We found that the GX 339-4 data 
can be fitted well without including any thermal component, and 
the entire broadband SED is well described by the superposition 
of synchrotron and inverse Compton photons from the pre-shock 
region near the jet base as well as optically thick synchrotrons 
from post-shock regions farther from the jet base. From the radio- 
IR slope, the required inclination is between 45-50 degrees, but 
not well constrained and hence was fixed at 47 degree. The length 
of the jet was fixed to 10'*^ cm and fraction of non-thermal elec- 
trons fixed to 0.75 for the same reasons given for XTE Jl 1 18+480. 
Since the spectrum of GX 339-4 could be entirely dominated 
by jet emission, an d the source was in a bright, hard X-ray state 
jHoman et al.l2005l) , the input jet power is much larger than that of 
XTE Jl 1 18+480. In fact the jet power for this obser vation is larger 
than any of the previous observations of this source jMarkoff et al.l 
I2OOII2OO5I) . Also the acceleration region starts much farther out 
for the GX 339-4 data set. 

If the temperature of the thermal electrons (T^) at the base is 



left completely free, as in "Model 1" (fifth column in Table|2ll, then 
Te for the best-fit model is comes out to be rather high, and sug- 
gests that the SED at energies > 10 keV is almost entirely due 
to SSC emission from the base (see top two panels in Fig. |4]l. 
The geometry of the base is rather compact and the source is in 
a bright hard state; therefore a strong SSC emission could imply 
that the photon density at the base becomes high enough for pair 
processes to become important. Since pair production is not in- 
cluded in the model, an estimate of the importance of pair pro- 
cesses for the best fit models were done as follows: the pair pro- 
duction rate (w^p) was calculated from the photon field at the base 
of the jet using the angle averaged pair production cros s section 
jGould & SchredeJI 19671 ; iBoettcher & Schli ckeiser"l997^ ; pair an- 
nihilation rate (hp„) was calculated following , S vensson (1982) . For 
"Model 1" of GX 339-4, hp,, ~ 10" cm-^j-' and hpa ~ 10'^ 
cm^'i"', suggesting that pair processes could be important in this 
case. Therefore we explored the parameter space, limiting the elec- 
tron temperature to be less than 5 x 10'" K and somewhat less com- 
pact base (larger ro, ho) to search for solutions which also satisfy 
the condition n,,,, < hp„ (see "Model 2", sixth column in Table |2] 
and bottom panels in Fig.|4ll. The fit derived number density (n), 
hpa and hpp for the different models presented in Table[2l are given 
in Table[l] 

In addition to the continuum jet model, an additional Gaus- 
sian line centered at 6.3 keV and reflection from the accretion disc 
is needed to obtain good fits for the GX 339-4 data. Presence of 
the line and reflection suggests the presence of an accretion disc. 
However, as in the case for XTE Jl 118+480, the disc is cold and 
therefore not detected by the PCA. 

The best-fit model parameters, associated confidence inter- 
vals, x~ statistic and corresponding chance probabilities for both 
sources are shown in Table |2] Given the simplicity of the model 
it is hardly unexpected that the reduced values obtained from 
the fitting are not very close to 1 in some cases. For instance, the 
viscous+irradiated disc model does not include details of disc at- 
mosphere or ionization balance equations, and is therefore not able 
to predict the spectral breaks in the optical data. Because of its sim- 
plicity, this is not meant to be a real model of the disc, although it 
incorporates the dominant ongoing physical processes and can re- 
produce the general shape of the broadband continuum. Moreover, 
while fitting a data set spanning 10 decades in frequency (or en- 
ergy) space, the long lever arm between radio and X-rays constrains 
the jet model more strongly than the overall chi-squared (see e.g. 
iMarkoff et ai] |2008'). Best fit models along with data and residuals 
are shown in Figures|2l [3]and|4l 



4 DISCUSSION AND CONCLUSIONS 

We have analyzed broad band multi-wavelength observations 
of the galactic black hole transient systems XTE Jl 118+480 
and GX 339-4, during their outbursts in 2000 and 2005 for 
XTE Jl 1 18+480, and during the 2002 outburst of GX 339-4. The 
conclusions from our study are summarized below. 

The data for XTE Jl 1 18+480 cannot be fit well with a jet con- 
tinuum model alone, due to an excess in the OIR fluxes. During 
these observations, in OIR the source was at least 3.5 magnitude 
brighter than q uiescence; t herefore the contribution from the Ml V 
secondary star jChatv et al. 2003, and references therein) is small. 
Modeling the data from the 2000 outburst and 2005 outburst of this 
source we found that the outer disc was strongly irradiated during 
the 2005 data. It appears that the inner edge of the accretion disc 
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was much closer to the central black hole during the observations of 
2005, thus providing an ample source of soft X-ray photons to irra- 
diate the outer disc. On the other hand, the sharp cutoff in the EUVE 
data and absence of any signature of irradiation from the outer disc 
suggests a recessed, cold accretion disc for the 2000 data. How- 
ever, as discussed in 33.3.21 the slope of the EUVE data depends 
very strongly on the assumed column depth, and a s mall decrease 
in Nh can flatten the EUVE SED significantly (see e.g. lHvnes et al.l 
120001) . Such a flat EUVE SED would b e expected from a disc with 
small rin and comparatively higher Tin jReis et alj2o"09h . 

Our calculations in 33.2.11 show that at any point on the sur- 
face of the accretion disc, the heating of the disc by the jet is seven 
or more orders smaller than the viscous heating. There are two fac- 
tors which largely reduce the amount of jet radiation that is incident 
on the disc: Firstly Doppler de-beaming - since the bulk motion of 
the jet is relativistic, special relativistic Doppler beaming concen- 
trates most of the emergent radiation from the jet in the direction 
of its motion, i.e. away from the disc. A second, but significant fac- 
tor that also contributes to the dimunition of jet flux incident on 
the disc is the fact that the post-shock jet synchrotron spectra peaks 
at increasingly longer wavelengths at increasingly larger distances 
from the base of the jet, thus creating fewer X-ray photons that can 
heat the disc. This is illustrated in Fig. [T] where we show the ra- 
tio of jet-heating to viscous heating as a function of radius. Since 
heating by the jet alone cannot produce the OIR excess in Fig.|2]we 
surmised that this the OIR excess is indeed caused by the more con- 
ventional form of irradiation, i.e. irradiation by soft X-ray photons 
from the iimer-disc. From an energetics point of view, it is easy to 
see that the inner accretion disc alone (see Fig.|2l at around 1 keV) 
radiates at least two orders or more energy than the outer region. 
Furthermore, it is well known that the outer edge of th e accretion 
disc can flare-up or warp (e.g. lOgilvie & DubuslboOlh . thus mak- 
ing the outer parts of the disc more easily visible from the inner 
regions. 

Fits to the 2005 data of XTE Jl 118+480 using the jet 
model+irradiated disc (described in 33.2.21 ) suggest a cold disc 
(rjn ~ 0.2 keV) with inner radius close to the innermost stable 
circular orbit (6 Rg for a Schwarzschild black hole, and smaller 
if spinning). Given the uncertainities associated with using a sim- 
ple T ~ R^ ^'"^ model which ignores co rrection factors like spectral 
hardening 

Ism 

mura & Takahara|[T99^ . full disc structure as well 
as pro per boundary conditions at the inner edge jZimmerman et al.] 
12005 ^. the errors associated with the estimates of the irmer disc 
parameters are most likely an underestimate. Another source of 
uncertainity comes from that fact that we assume the jet to be 
perpendicular to the disc. A misaligned jet jMaccarori3 1 20021) . as 
seems likely for XTE Jl 118-1-480, could introduce additional er- 
ror in the determination of r,„. Along with the uncertainities in 
instrument calibration, modeling and jet-disc alignment, the sharp 
low energy cutoff of PCA sensitivity below 3 keV makes detection 
of a cold disc with small ri„ even more difficult than a recessed 
cold disc. Nevertheless the numbers for r^^ and Ti„ we obtained 
for XTE Jl 118+480 are similar to those obtained by several oth- 
ers for black hole transie nt s in their hard state (e.g.loi Salvo etal 
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2007l : iReis et alj|2009l) . If such a cold disc with small ri„ is indeed 
realised, then Compton scaterring of the disc photons could be 
strongly anisotro pic, with most of the sc attered emission beamed 
toward the disc jHenri & Petruccilll997h . This could compensate 
partly the beaming effects due to the relativistic motion of the gas, 
and increase the illumination of the disc (both irradiation and re- 
flection effects). Also in this case there could be a radiative feed- 



back due to the local illumination of the disc by the base of the jet 
which i s similar to the situation i n standard accretion disc corona 
models ( iHaardt & Maraschill 1993 1) which would act as a thermostat 
keeping the temperature at the base lower than what is assumed 
here. A somewhat more conservative model, where the tempera- 
ture and radius of the inner edge of the accretion disc were fixed to 
0.1 keV and 30 Rg respectively (Model 2 in Table[2l( also describes 
the data well, showing the difficulty in detecting a cold disc with 
small ri„ from the available data. 

The best-fit outer radii of XTE Jl 118+480 for both Model 
1 and Model 2 are smaller than the circularization radiu s, using 
the la test values of orbital parameters published by Geli no et aU 
( l200d) . These results are therefore consistent with the assumption 
that accretion is occuring via Roche lobe overflow. The observation 
was made during the decline from outburst peak. Hence the dis- 
crepancy between the circularization radius and Ro,,, inferred from 
fits, in the context of standard disc instability models, could be due 
to an inward moving cooling wave. A decr easing trend in the o uter 
radiative disc radius was also observed bv lHvnes et all ( l2002h for 
the system XTE J1859+226. However the value of R„,„ in our case 
depends largely on the five OIR data points, all of which lie on the 
Rayleigh- Jeans region of the irradiated disc. Given the plethora of 
spectral f eatures seen in OIR spectra of X-r ay binaries in outburst 
(see, e.g.. lHvnesll200^ ; ICharles & Coell200Q) . the limitation of us- 
ing a few bands is obvious. 

Analysis of t he data fr om the 2 000 outbu rst of 

XTE Jl 118+480 by lEsin et al.l ( 1200 ih and lYuan et al.1 j2005l) 
showed that the UV and higher energy data can be explained by 
an ADAF model. The flux from the ADAF model alone however 
falls sharply at energies below optical and underestimates IR and 
any lower e nergy data (e.g. ra d io), thus requiring an additio nal jet 
component i lYuan et alj |200^; iNaravan & McClinto"c3 I2OOS) . In 
the iet-ADAF model iYuan et alj|2005h r the X-ray photons have 
a thermal Comptonization origin, whereas in the current work, 
the X-ray photons can originate from either synchrotron (e.g in 
XTE Jl 18+480), or inverse Compton scattering, or a combination 
of both (e.g. in GX 339-4). Better simultaneous broadband 
coverage and future X-ray polarization studies can help to estimate 
the relative importance of these two emission components which 
dominate the X-ray spectra of X-ray binaries in hard state. 

From the available data, the X-ray spectrum of 
XTE Jl 118+480 during both outbursts show little or no cur- 
vature at energies higher than ~ 10 keV. Reflection features like 
the iron line complex and the reflection "bump" near 20 keV are 
also extremely weak. Detailed analysis of the X-ray spectra taken 
by Chandra and RXT E, during the 2000 outburst of this source 
bv lMilleretZH2002h also showed extremely weak reflection fea- 
tures, as would be expected for a jet+recessed disc scenario. From 
a phenomenological point of view, since the curvature is small, 
using a broken power law to describe data of X-ray binaries in hard 
sta te has been shown to work remarkably wel l, e.g., for GX 339- 4 
bv 'Nowak et alj j2005l ) and for Cyg X-I bv lWilms et al] ( |2006|) . 
Residuals of a power law fit to the X-ray continuum in the 2005 
data set of XTE Jl 1 18+480 show small systematic variations near 
5-7 keV, which improves slightly upon addition of a Gaussian 
line. However, given the weakness of the line, its energy or width 
cannot be constrained well. Total contribution from various source 
of inverse Comptonization (SSC within the jet base and external 
Comptonization of soft thermal photons from the accretion disc) 
is also much smaller compared to that of GX 339-4. The X-ray 
region of the SED of XTE Jl 1 18+480 is dominated by post-shock 
synchrotron photons for both 2000 and 2005 outbursts. We note 
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that the electron acceleration region seems to be in the nozzle 
of the jet in the case of 2005 outburst data of XTE Jl 118+480 
while it is farther out in the jet in GX 339-4 as well as the 2000 
outburst of XTE Jl 118+480. While exhaustive searching of the 
space confirms that this is required by the data, it is not clear 
why the acceleration region is inside the nozzle for the 2000 data 
of XTE Jl 118+480. Another unexpected result is that the base 
of the jet seems to be quite small (ho/ro ~ 0.4) for the 2000 
outburst data of XTE Jl 118+480, whereas ho/ro ~ 1.5 for other 
outbursts of XTE Jl 1 18+480 as well a s other sources ( see e. g. 
iMarkofTetalJ ( |2005|) ; lOallo et al.l ilOOt : iMigliari et al.l tOOi) ). 
Given the uncertainties in our understanding of jet formation and 
the acceleration region, and the simplicity of the model, here we 
report the best fit parameters as obtained. The particle energy 
distribution index (p) of the relativistic post-shock electrons for 
data from both outbursts of XTE Jl 118+480 is ~ 2.5-2.6. This 
value is much steeper than what could be expected from an initial 
distributio n of particles accelerate d by a diffusive shock process 
(~ 1.5 - 2: lHeavens & Drurvll988h . However it is possible that for 
XTE Jl 1 18+480 the observed X-ray bandpass is at a higher energy 
than the cooling break energy, where the particle distribution 
steepens from p to p + 1 due to enhanced cooling. One of the 
key assumptions in our model is time independence of particle 
distributions. However if the steep power law drop-off of the 
observed photon index is indeed due to enhanced cooling, or more 
generally a departure from the assumed equillibrium between the 
heating and cooling processes, a more thorough time dependent 
analysis of particle evolution is neccessary. We will present time 
dependent particle evolution in the context of jets in a forthcoming 
paper. 

The jet inclination of XTE Jl 118+480 is still a matter 
of c oncern. The optically measured system inclination is near 
70° jGelino et al.ll2006l ; IChatv et alj|2003l : iMcChntock 61^1200 ll : 
IZurita et aL 2002). However, given the flat radio-to-IR fluxes, 
jet inclinations higher than ~ 30° would be very difficult to 
fit which might sugg est a large misalig nment between between 
the disc and the jet JMaccarond 12002 ). Evidence of misalign- 
ment between disc and jet has been seen i n several stellar X- 
ray binary systems, e.g., GRO J1655-40 tareene et al.' 200 ll: 
Hjellming & Rupen, 1 19951) , SAX J1819-2525 (Hiellming et al, 
XTE J1550-564 faannikainen et al 
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as well as in s everal AGNs, e.g., N GC 
3079 (Kondratko et al. 200 l), N GC 1068 jCaproni et al.l2006l) and 
NGC 4258 aCaproni et al]|2007 h. Allowing the inclination to vary 
for XTE Jl 118+480 during the fits gives a best-fit inclination of 
~ 25°, and it is almost impossible to find statistically good fits (e.g. 
with;t^^/v < 2) for inclinations larger than 30°. Previous fits by 
iMarkoff et all ( 1200 1[) also required a similar, small inclination. In 
our model we assume the lateral expansion and longitudinal accel- 
eration of the jet to be driven solely via the solutions of the relativis- 
tic Euler equation for adiabatic expansion. If this is not the case, and 
e.g. the jet radius and acceleration are determined by some other 
process(es) leading to a stronger magnetic field and/or smaller bulk 
velocity along the jet, that could lead to a fl atter radio-IR spectr um 
as well (also see the discussion in §3.3.3 o f lMigUari et ai]|2007l) . 

The data for GX 339-4 does not formally require a thermal 
component (e.g. an accretion disc or emission from the secondary 
star) and the entire broad band emission from radio to X-rays is 
satisfactorily well fit by non-thermal photons originating from pre- 
and post-shock synchrotron, and synchrotron self Compton emis- 
sion m ainly from the j et bas e. Previous fits to the 3-100 keV X-ray 
data bv lHoman et alj j2005l) also suggested that the disc contribu- 



tion within this energy range is ~ 1% of the total luminosity or less. 
However presence of the line near 6.5 keV and signature of reflec- 
tion suggests a cold accretion disc near the black hole which cuts off 
exponentially at an energy significantly below the lower detection 
limit of the RXTE/PCA. Since the X-ray spectrum has a luminosity 
of ~ 5% LEdd> the disc luminosity is presumably ~ 0.05% Leim or 
smaller. 

The data for GX 33 9-4 were tak en when the source was in a 
bright. X-ray hard state jHoman et al.u2005). The jet p ower Nj is 
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much higher than that obtained by|] 
led data from previous outbursts of the same source. We also obtain 
a higher equipartition factor (k) and shallower particle distribution 
index. It is interesting to note that ro and Zacc for GX 339-4 are 
always found to be higher than those found for other sources. If the 
electron temperature at the base is left unconstrained (e.g. "Model 
1" in the fifth column of Table |2j, then the density of photon field 
from best-fit model is high enough that pair processes may become 
important. In this case many photons are up-scattered above the 
pair production threshold, and the rate at which pairs are produced 
is much larger than the rate at which they are annihilated, thereby 
changing the density at the base of the jet. Constraining the elec- 
tron temperature to be less than 5 x 10'" K, we present another fit 
("Model 2", sixth column of Table |2]( where the photon density at 
the base is small and pair production rate is much smaller than the 
pair annihilation rate. For this fit the the fractional pair annihilation 
rate («pa/n) is also much smaller than unity, so that pair processes 
are not important for this fit. 

Using the irradiated disc+jet model presented in this paper 
we have shown that in certain cases (e.g. the 2000 outburst data 
of XTE Jl 1 18+480) there are not enough X-ray photons from the 
inner regions to irradiate the disc appreciably. The 2005 outburst 
data of XTE Jl 118+480 presents an intermediate case where the 
contribution from an irradiated outer disc is comparable to the jet 
emission in OIR. The 2002 data of GX339-4 presents the extreme 
case where the entire broadband SED is dominated by emission 
from the jet. 

The jet model has been used to analyze broad band SEDs of 
the stellar mass black h oles XTE 111 18+480 jMarkoff et iJ]|200ll 
this work), G X 339- 4 (Markofl' et al. 2005, this work), Cyg X-1 
tMarkoff' et al]|2005l). A 0620-00 (Galloetal. 2007J, GRO J1655- 
40 dMigUari et alJI 20071) . as well as that of th e supermassive black 
holes in the nucle us our own galaxy (S gr A*: lMarkoff et al.ll2007h 
and that of M81 ( iMarkoff et al]|2008h . From the analysis of the 
spectra of these sources we note that the geometry of the jet base, 
scaled in units of gravitational radius, does not vary largely. The 
radius of the jet base, Rq, lies in the range of 10-100 Rg, and the 
ratio of ho/ro at the nozzle lies within 0.4—11, both suggesting a 
relatively compact jet base and the universal nature of mass scal- 
ing in accretion onto compact objects. The fraction of particles ac- 
celerated to a power law distribution beyond Zacc, while not very 
well constrained by the fits, preferes a value in the range of 0.6-0.9 
and is fixed to 0.75. Considering that most of the known accelera- 
tion mechanisms in astrophysics do not have such a high efficiency, 
this number seems rather high and requires further exploration. The 
particle acceleration index (p) lies within 2.2-2. 7, which is steeper 
than that expected from a relativistic shock ^Heavens & Drurvl 
Il988l 1.5-2, see, e.g.), suggesting that in these cases the X-ray en- 
ergies may lie beyond the cooling break ( Kardashev 1962), so that 
the power law index of observed SED steepens by an additional 
amount of 0.5. 

The magnetic-to-particle energy densities at the jet base (k) is 
often found to be greater than unity. MHD simulations of relativis- 
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tic jet launching i McKinnev 8i GammijI 20041: Nakamura & Meie3 
|2004 iDe Villiers et alj|2005l : iFragile & Meieill2008h suggest that 
near the jet base, close to the launching points, the jets are mostly 
Poynting flux dominated and have small plasma fractions (ratio of 
gas pressure to magnetic pressure is less than unity). We are fo- 
cusing on inner regions where this is likely still to be the case. In 
GRBs and AGN, most data is only relevant to regions well beyond 
the magnetosonic fast points, where the flow has been driven to- 
wards equipartition via conversion of magnetic energy to particle 
energy. 

Temperature of the thermal electrons at the jet base (Tg) is in 
the range of 2 - 7 x 10'" K for the stellar black holes and ~ 10" K 
for the supermassive black holes. The flow near the jet base could 
be a radiatively ineflicient accre t ion flow (RIAF), i. e. an ADAF 
jNaravanetaiJl 19951 : lReej|l982l : Ishapiro et al.lll976h, o r in light 
of the the magnetic domination, an MDAF ( Meiei 20051) . In such 
a case, the ion temperature ca n be as high as the virial tempera- 
ture, and electrons at 10'"'" K dShapiro et alJI 19761 : iNaravan et al.l 
1 19961) . Keep in mind that for an outflow model, the particles are 
also not required to be gravitationally bound, and processes in the 
jet launching can heat the particles beyond virial, in theory. 

Most likely the similarity in the derived physical parameters 
among the diverse sources and their diverse accretion states is not 
coincindental, and may hint towards an inherent similarity in the 
process of formation and propagation of compact jets. The results 
presented above show the importance of quasi-simultaneous, multi- 
wavelength, broad band SEDs in constraining geometrical and ra- 
diative parameters associated with accertion flows near compact 
objects. 
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XTE Jl 1 18+480 (2000 outburst) 


XTE Jl 1 18+480 (2005 outburst) 


GX 339-4 (2002 outburst) 


MJD of observation 


51652 


53393 


52367-68 


Spectral Coverage 








Radio: 


15 GHz 


15 GHz 


5 GHz 


UV,Opl/IR: 


M, L, K, H, J-bands, EUVE and HST spectra 


K, H, J, V, B-band 


H, I, V-band 


X-ray: 


Chandra (0.24-7 keV), RXTE (3-1 10 keV) 


RXTE (3-70 keV) 


RXTE (3-200 keV) 


A'h (10^2 cm-^) 


0.013 


0.013 


0.6 


Mbh (Mo) 


8.5 


8.5 


7 


Distance (kpc) 


1.7 


1.7 


6 



Table 1. Spectral coverage and source parameters used for the fits. For the 2000 outburst of XTE Jl 118+480 the radio data were obtained using the Ryle 
telescope, IR using the UKIRT, optical spectra u sing the HST, UV spectra using EUVE, X-ray spectra using Chandra and RXTE. The SED was constructed 
from data published in lMcClintock etaljj200lh . For the 2005 outburst of XTE Jl 118+480 the radio data were obtained using the Ryle telescope, IR using 
the UKIRT and opt ical using the Liverpool telescope. The column density, mass of the black hole and distances we re taken from the rec ent observations by 
iGelino et all j2006l) for the source XTE J l 1 1 8+480. For GX 339-4, the SED was constructed fror n data pubhshed inlHoman et alj )2005l) , black hole mass is 
lower limit bv lMunoz-Darias et alj j2008l) , distance and column density from the ranges given by faynes et alj2004h . 
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Parameter (Unit) 


XTE Jill 8+480 

(2005 outburst; Model 1) 


XTE Jill 8+480 

(2005 outburst; Model 2) 


XTE Jill 8+480 
(2000 outburst) 


GX 339-4 

(2002 outburst; Model 1) 


GX 339-4 

(2002 outburst; Model 2) 


N, (10-3 LEdd) 


9 ^9+0.02 
-0.06 


9 77+0.01 
-0.05 


7 9+0.1 
'■ -0.5 






'-0 (Rg) 


a on+0.18 
'^■-^^-0.17 


11.8!«i 


20.8!]° 


33!3 


114+], 




1 CQQ+O.OB 

-^°'-0.002 


1 ^n+0.03 
-^'-0.02 


0.4* 


1 7+"-2 
^■'-0.1 


1.51+°;^^ 


Te (10'° K) 


A. 99+0.11 

-0.02 


4 04+" °' 
^■"^-0.02 


3.76+J];g? 


6 7+"" 

"■'-0.5 


^■^-02 


Zacc (Rg) 


12 6+'^ 


12 6+2 


8.2+» 


252+212 

^-^ -130 


526!«5 


P 


2 63+" '" 

^■"^-0.01 


9 71+0.01 
^■'^-0.01 


2 48+" 02 
"^■^°-o.oi 


9 9Q+0.03 


2 20+" "2 

^■^"-0.01 


Esc (10-3) 


Q+0.02 
'■'-0.1 


90 Q+0 3 




16+""'* 
"■'"-0.08 


44+2 


k 


I 9+0.1 


2 3+°-' 


2 1+°' 


9 9+1.4 
^■^-0.8 


9 qo+0.01 
2.iS„m 


Inclination (deg) 




30* 


30* 


47* 


47* 


nn (Rg) 




30* 


341+^2 


Not constrained 


Not constrained 


Tin (keV) 


9^4+0.004 


0.1* 


025+" ""' 


Not constrained 


Not constrained 


'out (Rg) 


15700+ 

-350 


8850+]™ 


29400;^.J^Q 


Not constrained 


Not constrained 


To,,, (K) 


18100!r 


3O5O0;jj;° 


650+-'^" 

-550 


Not constrained 


Not constrained 


/(lO-'t) 


1.2 


1.0 


Not constrained 


Not used 


Not used 


^Lmc (keV) 


6.63* 


6.63* 


Not used 


6 3+"-2 

"■^-0.2 


6 3+"-3 

"■^-0.2 


WLine (eV) 


39.3 


27.7 


Not used 


93.5 


58.7 


a/In 


Not used 


Not used 


Not used 


22+" '" 

"■^^-0.14 


35+" "2 

■-^ -0.02 


x^/y [Q] 


70/50 [0.032] 


90/53 [0.0011] 


140/137 [0.406] 


137/112 [0.054] 


157/112 [0.0032] 



Table 2. Best fit parameters for XTE Jl 1 18+480 and GX 339-4, and Ax^ < 2.71 confidence intervals. A star (*) next to a number indicates that the paiameter 
was frozen to that value. Nj = jet normalization; ro = nozzle radius; fto/ro = height-to-radius ratio at base; = pre-shocic electron temperature; Zacc = location 
of acceleration region along the jet; p = spectral index of post-shock electrons; e^c = 0.36/ratio of scattering mean free path to gyroradius (see 33. It ; k = ratio 
between magnetic and electron energy densities; Inclination = orbital inclination as seen from Earth (also the jet inclination); rjn/out = radius of the inner/outer 
edge of the irradiated accretion disc; Ti^/oat = temperature at the inner/outer edge of the irradiated accretion disc. While Ti,„r is in Kelvin, T[„ is reported in keV 
as is customary for X-ray observers; / = (2;r)-'x the sohd angle subtended by the outer, in'adiation dominated disc as seen from the inner disc; Eunc = energy 
of the iron Kcv line complex from the gaussian model; WLinc = equivalent width of the line; Q./2n = reflection fraction from the reflect model, x^/^ is the 
reduced chi-squared and Q is the coiTespond ing chance probability. For fitting the data of XTE Jl 1 18+480 from 2000, the energy range of 0.15 - 2.5 keV was 
excluded (see text, and also lEsin et alj200lh . Also the disc inner edge parameters (rj^, Tj^) are extremely sensitive to the extinction column depth ( j|3.3.2t . 



Source/Outburst/Model 


log (fi [cm 3] ) log 


{hpc, [cm 3 s ' 


] ) log (hpp [cm 3 s '] ) 


XTE 111 18+480/2005/1 


13.9 


12.0 


7.6 


XTE 111 18+480/2005/2 


13.8 


11.8 


6.8 


XTE 111 18+480/2000 


13.7 


11.7 


5.8 


GX 339-4/2002/1 


14.2 


12.4 


17.1 


GX 339-4/2002/2 


13.4 


11.0 


6.1 



Table 3. The lepton number density (n), pair annihilation rate (hpa), and paii' production rate (hpp) at the base of the jet for the various models presented in 
Table[2] 



12 D. Maitra et al. 



I ' — ' — I ' — ' — I ' — ' — I 




J I I ■ ■ ■ I I I I I ■ ■ ■ ■ I I I I I ■ ■ ■ ■ I I I I I I ■ ■ ■ I 

10^ 10^ 10^ 10"^ 



Figure 1. Comparing viscous heating (//vise) ^ind maximum possible "jet-lieating" (//jet) for XTE Jl 1 18+480. The filled circles are the computed upper limits 
of the ratio of jet induced heating to viscous heating (//jet///visc), as discussed in ^3.2.11 joined by a smooth spline. /"]„, Rjn and /?„„, for the disc were 
determined from broadband "Model 1" fits to the SED. For R > 100 Rg, the jet induced irradiation heating term falls more slowly (~ R-- *) than the viscous 
heating (~ /?"'), which causes the ratio //jct///visc to slowly increase at increasing radii. For fiin < ro, the jet heating saturates but the viscous heating continues 
to increase for smaller radii, causing //jci///visc to drop sharply for decreasing radii. 
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Figure 2. Jet model fits and residuals for the XTE Jl 1 18+480 data from its 2005 outburst. Left panels show the entire broad band SED from radio through 
X-rays, and a zoom of the X-ray region is shown in the right panels. Radio and IR data from the Ryle telescope and UKIRT respectively are shown by purple 
squares, RXTE/PCA data are shown by by blue circles and RXTE/HEXTE data by orange triangles. The top panels are for "Model 1" in Table |2] (second 
column) where along with other parameters (see text), the jet inclination, ri„ and Tin, were also allowed to vary. The value of the reduced chi-squared for 
best fit parameters of "Model 1" is 1.41. The bottom panels are for "Model 2" in Table |2] (third column) where jet inclination=30° , rin=30 Rg and Tin =0.1 
keV. The value of the reduced chi-squared for best fit parameters of "Model 2" is 1.69. The dark-green dash-dotted curve shows the post-shock synchrotron 
contribution, the green dashed curve shows the pre-shock synchrotron, blue dotted curve shows the Compton upscattered component. Flux from the irradiated 
accretion disc is shown by the orange dash-dotted curve. The solid grey hne shows the total jet+disc model continuum spectrum without convolving through 
detector responses, interstellar extinction, the iron line or reflection. The iron line near 6.6 keV is shown by the thick, solid green line. The red line shows 
the properly forward-folded models taking into account detector responses, interstellar extinction, the iron line and reflection. Since interstellar extinction is 
small for XTE Jll 18+480, and the iron line and reflection features are very weak too, the forward-folded model (red line) is almost indistinguishable from the 
jet+disc model continuum (solid grey) in this figure as well as in Fig. [3] However for GX 339-4 (Fig.|4), where the iron line as well as reflection features are 
stronger, and the extinction is also larger, the dift'erence between the unfolded and folded models (grey and red lines) becomes apparent. 
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10"" 10"^ 10"^ 10"^ 10"* 10"-' 0.01 0.1 1 10 100 2 5 10 20 50 100 



Energy (keV) Energy (keV) 

Figure 3. Jet model fits and residuals for the XTE J1118+480 data from its 2000 outburst (taken from' McClintock et al]l200lh are shown in the left panel. 
The radio data from the Ryle telescope and IR data from UKIRT ai'e shown by red circles, HST spectrum using green, EUVE data using orange triangles, 
Chandra spectr um using blue ci rcles and RXTE data using purple squares. Note the dip in the 0.15-2.5 keV region which is attributed to the presence of a 
warm absorber ( Esin et alj200ll) . This region was excluded from our fits. A zoom of the X-ray region after excluding 0.15-2.5 keV data is shown in the right 
panel. The colour coding for the model components are the same as in Fig.|3] 
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10"° 10"' 10"° 10"' 10"* 10"' 0.01 0.1 1 10 100 2 5 10 20 50 100 



Energy (keV) Energy (keV) 

Figure 4. Same as Fig. [2] but for the source GX 339-4. The models presented here are purely jet models without any contribution from an accretion disc. The 
top panels are for "Model 1" (Table|2] fifth column) where along with other parameters the temperature of the thermal electrons at the base of the jet (T^) was 
allowed to vary freely. The reduced chi-squared for this model is 1.22. Note the relatively strong inverse Comptonization component in the X-rays, originating 
at the jet base due to SSC emission in this model, compared to that in the case of XTE Jl 118+480. The bottom panels are for "Model 2" (Table|2l sixth 
column) where we constrained Te to be less than 5 X lO'" K, reducing the SSC emission at the base. The base is less compact and acceleration starts farther 
out along the jet for "Model 2" compared to "Model 1". The reduced chi-squared for this second model is 1.40. While the first model is a statistically better fit 
than the second, the photon density at the base can become high enough for pair processes to become important in the first model (see 33.3.3| j|4land Table[3), 
which is not the case for the second model. Also note the strong iron line near 6.3 keV and the reflection "bump" near ~ 20 keV, suggesting the presence of a 
cold accretion disc not detected at energies to which RXTE is sensitive. 
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